Peak Effect in Superconductors: Absence of Phase Transition and Possibility of 

Jamming in Vortex Matter* 
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(Dated: February 2, 2008) 

The magnetic field _B dependence of the critical current / c for the vortex phase of a disordered 
superconductor is studied numerically at zero temperature. The I C (B) increases rapidly near the up- 
per critical field B C 2 similar to the peak effect (PE) phenomenon observed in many superconductors. 
The real space configuration across the PE changes continuously from a partially ordered domain 
(polycrystalline) state into an amorphous state. The topological defect density nd{B) ~ e aB with 
k > 1 for B > 0AB C 2- There is no evidence of a phase transition in the vicinity of the PE suggesting 
that an order-disorder transition is not essential for the occurrence of the PE phenomenon. An 
alternative view is presented wherein the vortex system with high dislocation density undergoes 
jamming at the onset of the PE. 

PACS numbers: 



I. INTRODUCTION 



The critical current density J c for a type-II supercon- 
ductor decreases monotonically as a function of the mag- 
netic field B or the temperature T over most part of 
the B-T phase diagram. But quite often, a peak in the 
J C (B) (or J c (Tj) is observed close to the upper critical 
field B c2 (T) (or critical temperature T C (B)). This phe- 
nomenon is known as the peak effect (PE) and was first 
studied in 1961. Extensive work since then have shown 
that PE exist in wide class of materials which includes 
low-T c (Nb and its alloysQ, V 3 Si0, NbSe 2 [3[), high-T c 
(YBa 2 Cu 3 7 |l), heavy fermion (UPt 3 @, UPd 2 Al 3 |3), 
and other exotic (Sr 2 Ru04[7]]) superconductors, includ- 
ing the recently discovered compound MgB 2 [j|. Apart 
from single crystals, the phenomenon have been observed 
in polycrystalline and thin amorphous films 0- The PE 
is found to be more pronounced in the weak pinning 
limit and progressively diminishes with increasing pin- 
ning strength[^J,[n|. An important characteristic of the 
PE phenomenon is the strong memory and history depen- 
dence of the various properties, and have been observed 
in static and dynamic experiments [p3 . Il3l ITU ITM lliil Il7| . 

The explanation for the PE is based on the softening 
of the vortex lattice with increasing B, first suggested 
by Pippard 18]. If the vortex lattice rigidity falls more 
rapidly than the elementary vortex pinning strength f p 
as B — > B C 2, the softening would allow vortices to con- 
form to the random pinning potential better compared 
to a stiff vortex lattice. This idea was further developed 
and quantified in the collective pinning theory by Larkin 



"This work was carried out during the authors' stay in the Depart- 
ment of Physics, University of California, Davis, 95616. 
t Present Affiliation: Materials Research Laboratory, John F Welch 
Technology Centre, GE India Technology Centre, Phase 2, Hoodi 
Village, Whitefield Road, Bangalore, 560066; Email address: ma- 
hesh. chandran@ge.com 



and Ovchinnikov|19j (LO). In the LO theory, the vor- 
tex lattice is broken into coherently pinned regions of 
volume V c = B? C L C within which the pinning forces act 
randomly. R c and L c are the transverse and the longitu- 
dinal correlation lengths, respectively, and are related to 
the tilt modulus C44 and the shear modulus Cqq of the 
vortex lattice. The balancing of elastic energy and the 
pinning energy within V c gives the critical current density 
Jc oc fpijip/Vc) 1 / 2 . LO suggested that in 3D the spatial 
dispersion of C44 can lead to an exponential decrease in 
V c near B C 2 giving rise to the PE phenomenon. The PE 
is not expected in the 2D since the only relevant elastic 
modulus 6*66 is weakly dispersive. 

Experimentally, the PE have been observed in thin 
amorphous films (samples with thickness d < L c ). This 
was interpreted within the collective pinning theory as 
due to dimensional crossover from 2D to 3D in the pin- 
ning characteristic [20j. Recently, the PE have been ob- 
served in ultra thin Bi films of d = 26 A°[2lJ which is 
difficult to account for within the dimensional crossover 
scenario. Moreover, neutron scattering experiment on Nb 
crystal shows that L c actually increases with B, and the 
decrease in L c in the PE region is not appreciable p^. 
Also, the LO theory cannot account for the PE in 
Jc(T)M- 

In the LO picture, the pinning energy overcomes the 
elastic energy beyond R c , and hence, the long range order 
is not expected to survive in the vortex lattice. In high- 
Tc superconductors, topologically ordered vortex lattice 
on a length scale r » R c was observed pil [gjj. De- 
tailed calculations showed that in the presence of weak 
point impurities a quasi-long range ordered vortex lat- 
tice survives in 3D with C(r) ~ r - ^, where 77 is a non- 
universal exponent ptj I27I ] . Such a phase has been termed 
as the Bragg glass (BG). The BG undergoes a first-order 
melting transition into a vortex liquid (VL) phase on in- 
creasing T. Also, the BG undergoes a weak first-order 
transition into the vortex glass (VG) phase on increasing 
j3pc|. In the VG phase, the plastic deformation caused 
by the topological defects leads to strong pinning, and 
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hence higher J c compared to the BG phase. The PE 
have thus been interpreted as a consequence of an order- 
disorder transition|29i l3Cj . Furthermore, the origin of 
anomalous dynamics near the PE is attributed to the 
coexistence of BG and VG phases in the vicinity of the 
transitionpul • 

Interpretation of the PE phenomenon as an order- 
disorder transition has limitations. Firstly, the BG phase 
is not expected in 2D|33|. and hence order-disorder tran- 
sition cannot account for the PE in amorphous films. 
Secondly, the BG phase have been shown theoretically 
to exist in the presence of weak point impurities. It is 
not clear if the BG would survive in polycrystalline sam- 
ples in which the PE has been well documented. Recent 
Bitter decoration of NbSe2 across the PE shows no dis- 
tinct ordered and disordered phases^- Also, the width 
of the field range in which the ordered and the disordered 
phases coexist have been shown to scale with the sample 
size, which is an evidence against the phase transition 
interpretation of the PEp^j. 

Attempts have been made to understand the PE phe- 
nomenon using numerical simulation. Cha and Fertig[3fij 
used an attractive term in the inter-vortex potential to 
tune the Cqq in 2D. A small peak in the total de-pinning 
force was observed which was associated with the tran- 
sition from the elastic to plastic depinning. But as em- 
phasized in Ref . [37l| . the peak in total depinning force as 
a function of Cqq does not imply a peak in J C (B). Later 
simulation in 3D showed an insignificantly small peak 
near the BG to VG transition 38] . The PE was also stud- 
ied in a model of layered superconductors where the peak 
in J c was induced through the decoupling transition j3fj. 

In this paper, the behaviour of the critical current 
I C {B) for a 2D disordered vortex matter is studied us- 
ing numerical simulation. The effect of finite vortex core 
is included in the inter-vortex interaction and is shown 
to be essential for the realistic description of the system 
close to B C 2- The I C (B) decreases monotonically in the 
intermediate field range and rises rapidly close to B C 2 
similar to the PE phenomenon observed in real systems. 
The vortex configuration shows no order-disorder tran- 
sition but a continuous transformation from a polycrys- 
talline domain state below the PE to a liquid-like amor- 
phous state in the PE region. An alternative scenario 
for the PE is discussed where the vortex system under- 
goes jamming due to increased dislocation density. An 
empirical relation between the I C (B) and the real space 
defect density is also obtained from the simulation which 
captures the essential physics governing the behaviour of 
I C (B) for a disordered vortex matter. 

The paper is organized as follows: in section 2, the 
simulation details are presented, followed by results in 
section 3. The results are discussed in Section 4, and the 
conclusions are summarized in Section 5. 



II. SIMULATION METHOD 

Consider a 2D cross-section of a bulk type-II super- 
conductor perpendicular to the magnetic field B = Bz. 
We model the vortices as classical particles interacting 
via two-body potential, and governed by an overdamped 
equation of motion 

rlr 

V^ = ~J2 VE7S(r< - r,) - ]T VE/'fc - R fc ) + F ext ,(l) 

j^i k 

where rj is the flux-flow viscosity. The first term on the 
left hand side represents the inter- vortex interaction, and 
the second term is the attractive interaction between the 
vortex and the quenched impurities. F ext is the Lorentz 
force on the vortex due to an applied current. Each of 
the three terms are further discussed in detafl below. 

The inter-vortex potential Uq(t) — s ^ 2 ° x2 K (f/X), 
where Kq is the modified Bessel function, and f = 
(r 2 + 2£ 2 ) 1 / 2 . 0o is the flux quantum, and A and £ are the 
penetration depth and the coherence length of the super- 
conductor respectively. This form of the interaction was 
derived by Clem from the Ginzbug-Landau (GL) equa- 
tion using variational method |40l |41| . Conventionally, 
the vortex dynamics have been studied using the poten- 
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tial U1{r) = 87 r2° A 2 Ko(r/X) which is derived from the 
London's equation. For small r, the K (r) oc — ln(r), 
leading to an unphysical singularity in U^(r) as r — > 0. 
The potential (r) is thus inadequate for describing cor- 
rectly the vortex core region which limits its applicability 
to small B. On the other hand, Clem's potential Uq(t) 
is well behaved as r — > and gives realistic form of the 
magnetic flux density B(r) and the supercurrent density 
J(r) near the vortex core region. The correct descrip- 
tion of the vortex core region is essential to extend the 
simulation to high magnetic fields where the cores of the 
neighboring vortices tend to overlap. 

An important consequence of using Clems's potential 
is that the vortices are not point particles but have a 
finite radius £. This is an improvement over the previ- 
ous simulations which treated vortices as point particles 
interacting via London's potential U"l(r). Note that the 
London's model of vortices interacting via two-body po- 
tential is strictly valid for B < Q.25B^-> |4l|. Even though 
Clem's potential overcomes some of the limitations of the 
London's potential at high fields, it still provides only a 
qualitative behaviour of the vortex system. A more re- 
alistic approach for understanding the vortex dynamics 
close to B C 2 is through the time dependent GL theory 
(TDGL). Unfortunately, the computation cost for TDGL 
is prohibitive even for a small sample size|42| which re- 
stricts its utility. 

The second term in Eq.(l) is the attractive force be- 
tween the vortex and the quenched impurities which 
act as pinning centers. The force is derived from the 
parabolic potential, U p (r) = Uo(^z — 1) for r < r p , 

and otherwise. The impurities are randomly placed 
at positions in the simulation box. The third term 
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is the Lorentz force F ext = x (foz where J is the 
the current density. A uniform current density along 
the y-direction is assumed such that all vortices expe- 
rience the same force F* xt along the x-direction. For 
the simulation purpose, the length is defined in units of 
X(B = 0) = Xq. The current density J and the velocity of 
the vortex v x are in the units of cfo/<t>a and fo/rj, respec- 

,2 

tively, where /o = gj^a ■ From the relation E = v x B, 
the v x oc E y oc V where V is the voltage generated in 
the direction of the current I oc J y = Ff xt . Hence, the 
v x (F^ xt ) behavior represents the V(I) characteristic (or 
equivalently E(J) curve) for a superconductor. 

The dimensionless magnetic field is defined as b = 
B/B C 2, where the upper critical field B C 2 = and 

£o = £(B = 0). The b is calculated from the lattice con- 
stant ^ = (i|)|(_L_)|. The GL parameter k = | is 
an input to the simulation. In the GL theory, the length 
scale A and £ increases with b and diverges as b — > 1 . The 
renormalization of A and £ with increasing b is included 
in the simulation through the relation X(b) — f(b)Xo and 
£0) = /(&)&, where f(b) = 1/^(1 - b 2 ). This form 
of f(b) is similar to the T-dependence of £ in the GL 
theory|43| with T/T c replaced by (B/B c2 ) 2 , and have 
been used previously in Ref.|4J|. The same form of f(b) 
is chosen for A and £ in order to keep k independent of 
B. 

The simulation was carried out by numerically inte- 
grating the equation of motion using the fourth order 
predictor-corrector scheme. Parallel algorithms were im- 
plemented details of which can be found in Ref.^f|. A 
perfect vortex lattice driven by a large current I 3> I c was 
used as the starting condition for the simulation. The / 
is then reduced to in small steps and the average volt- 
age V was calculated in the steady state at each step to 
obtain the V(I) curve. In the dimensionless units defined 
above, the V = I in the asymptotic limit / — > oo. The I c 
is defined as the current at which the V < 10~ 5 ^(J. The 
procedure used for obtaining the V(I) curve was moti- 
vated by experiments where a stable vortex configuration 
is observed if the system is brought to rest after driven 
with / 3> I c 0|- Also, the I c (b) from such a method 
shows the largest peak compared to that obtained from 
the field-cooling experiments. 

The real space configuration of vortices was charac- 
terized by locating the topological defects in the system 
using Delaunay triangulation. For the triangular vortex 
lattice, a topological defect is a vortex with coordination 
number other than 6. The defect density n^/Ag is defined 
as the number of defects per unit area of the simulation 
box. Most of the defects form dislocations which are 
bound pair of disclinations (vortices with coordination 5 
and 7). The fraction of free disclination is small over a 
wide range of fields, and hence rid ~ ^ridisi where ridisi is 
the dislocation density. Along with the V(I) curve, the 
behaviour of rid(I) was also determined for each value of 
b. The rid discussed below represents the defect density 
at I = 0, unless specified otherwise. 



The parameters used in the simulation are k = 10 and 
X(B = 0) = 1000A These values are typical of low- 
T c superconductors, particularly NbSe2- The periodic 
boundary conditions were imposed in both directions. 
The magnetic field b was varied by changing the size of 
the simulation box, keeping the number of vortices Af,, 
fixed (for smaller system size, N v is also allowed to vary 
between 800-1200). The results presented below are for 
N v = 4096. The prefactor Uq of the pinning potential 
U p (r) is distributed randomly between A ± 0.01 where 
A = (Uq). The simulation is restricted to point impu- 
rities with the range of the pinning potential r p = £o- 
The pin density n p = 2.315/Aq. The dilute limit of n p is 
chosen keeping in mind that the PE is observed close to 
B C 2 where the concentration of effective pinning centre is 
expected to be very small. 
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FIG. 1: The critical current I c (b) for A = 0.02 and N v = 
4096. The I c (b) for A = 0.04 is shown in the inset. For this 
plot, the N v was also varied between 800 and 1200 along with 
the simulation box size. The thick line represents the fit to 
the simulation data (shown by open symbols). 



III. RESULTS 
A. Critical current I c (b) 

Figure 1 shows the I c (b) for A = 0.02. For b < b p . on = 
0.8, the I c decreases monotonically. At b p , n, the I c (b) 
turns around and rises rapidly till b = 0.925. The rapid 
increases in I c above b pxm signals the appearance of the 
PE. The inset shows the PE for A = 0.04 for which a 
smaller system size was used. The b p , on decreases from 
a value of 0.9 for A = 0.01 to 0.75 for A = 0.04. The 
proximity of the PE to the upper critical field B c2 is in 
good agreement with the experiments. 

For b > b p = 0.925, the I c decreases. For b > 0.95 
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two vortices occasionally come together to form an ef- 
fective vortex with charge 2(f> . Overlapping of the two 
vortex cores is strongly influenced by the impurities. This 
was confirmed by simulating without the pinning centres 
where no overlapping of the vortices was found up to 
b = 0.98. The effective interaction between a 2<f>o vortex 
and the <po vortex is greater than the interaction between 
4>o vortices. This enhances the local rigidity of the vor- 
tex system, and possibly decrease the I c . It is worth 
mentioning that a short range attractive interaction ex- 
ist between the vortex cores 41] which is not included in 
the simulation. The important point is that the overlap- 
ping of vortex cores is observed for b > 0.95, and hence 
does not influence the increasing branch of the PE. It is 
possible that the I c (b) above b p is influenced by the over- 
lapping of vortices, a issue which can be addressed only 
by solving TDGL equations in the presence of impurities. 

Attempts were made to fit the I c (b) data to a single 
function across the PE. A reasonably good fit was ob- 
tained using the expression 

J c (b)=J c0 exp(-/36 fe ) 1 (2) 

with k w 2.6 and (3 w 9.75 for A = 0.02. For A = 
0.04, k « 2.4 and (3 w 8.3. The above function was 
motivated by the behaviour of the nd(b) which follows 
the same exponential form, whereas the dependence on 
(i\a<) reflects the behaviour of A and £ as b — ► 1. The 
renormalization of the length scales A and £ leads to weak 
inter-vortex interaction, and explains why the PE in real 
system occurs close to the B C 2(T) line. As shown below, 
the bare length scale Ao and £o leads to a stiff vortex 
lattice for b > 0.6 and consequently the PE is absent. 



B. V(I) characteristics across the Peak Effect 

Figure 2 shows the V(I) curves across the PE. The 
curves for high-6 crosses the curves for \ow-b in the sub- 
ohmic region. The crossing can be understood from the 
dynamics close to I c . For I > I c , the dynamics is hetero- 
geneous with few vortices moving in channels in a back- 
ground of defective vortex configuration^)^,!^. The 
width of the channels ~ ao and display large transverse 
wandering relative to the direction of the Lorentz force. 
Since not all vortices are active, the voltage response is 
sub-ohmic. The number of active channels increases with 
I, and for I — I p all vortices move in channels which 
are ordered transverse to the flow direction. The time- 
averaged defect density shows an abrupt drop at I p , and 
represents annealing of the plastic deformation (some dis- 
locations with Burger vector parallel to the Lorentz force 
can exist for / > I p ). For I > I p , the response is ohmic 
and the resistance approaches the free flux-flow value. 
The sub-ohmic response window 5Idi S = Ip — Ic thus 
quantify the extent of plastic deformation of the vortex 
system. As shown in Ref.|5lj|. 5Idi S (b) increases rapidly 




FIG. 2: The V(I) curve across the PE. Arrows point to the 
crossing of V(I) curves. Inset: the V(I /I c (b)) curves for (a) 
b < b P:0n , and for (b) b > b P:0n . The arrows represent the 
direction of increasing b. The magnetic fields in (a) are 0.3, 
0.4, 0.5, 0.6, 0.65, 0.7 and 0.75. In (b), the 6 increases from 
0.8 to 0.9 in steps of 0.025. 

in the PE region implying that larger currents are re- 
quired to anneal the defects. This is reflected in the slow 
growth of V(I) curve for high fields, which consequently 
cross the curves for low fields as shown in Fig. 2. 

The inset in Fig. 2 shows the change in V(I) charac- 
teristics as b is varied across the PE. For each curve, the 
I is scaled by the respective I c (b). Far below the PE, 
the curves are convex for I > I c . As 6 approaches b p>on , 
the curvature changes to concave. The concave part of 
the V(I) curve indicates a slow growth of the response 
due to increased plasticity of the system. Interestingly, 
for I/I c (b) > 1, the V(b) reflects the non-monotonic be- 
haviour of I c (b). Thus, the PE is characterized by static 
as well as dynamic changes in the system, and agrees 
with the conclusion from the experiments 52] . 



C. Real space configuration 

The simulation allows the changes in the real space 
configuration to be characterized precisely, thus allowing 
theories based on configurational changes to be tested. 
Fig. 3 shows the real space images at I — in a re- 
gion of the simulation box as b is changed across b ViOTL . 
The defects are marked by filled circles. A topologically 
ordered vortex lattice is not observed for any value of b 
consistent with the theory for 2D system 26] . But par- 
tial order can still be seen in large domains for b between 
0.3 — 0.6. The dislocations are arranged in string-like 
structures and forms the grain boundary (domain walls) . 
There are no free dislocations inside the domain, and free 
disclinations are absent below b PiOTl . The configuration in 
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FIG. 3: The Delaunay triangulation of the real space con- 
figuration across b v , on . The &=0.70 (a), 0.75 (b), 0.80 (c), 
0.85 (d), 0.875 (d), and 0.90 (f). The red and the blue dots 
represents vortices with 7 and 5 neighbors, respectively. The 
vortices with 4 and 8 neighbors are denoted by the green dot. 
The pinning strength A = 0.02 and N v = 4096. Only a small 
region of the simulation box is shown for clarity. 



the intermediate fields is reminiscent of a polycrystalline 
solid. With increasing 6, the domain size decreases but 
the increase in the defect density occurs within the grain 
boundaries. Some of the domains can also be seen for 
b = 0.70 and 0.75. The average size of the domains 
decreases from w 6ao below the PE to w 3ao on the 
increasing branch of I c {b). At b p , the configuration is 
amorphous (Fig. 3(d)) with defects forming a dense ho- 
mogeneous network. It is not possible to isolate individ- 
ual dislocations and the configuration is that of a frozen 
liquid 54] . 

The continuous transformation of a less disordered 
state below 6 Pj0ra into an amorphous state at b p pre- 
cludes an order-disorder (or BG to VG) transition un- 
derlying the PE. The PE is thought to be due to such 
a phase transition in the vortex system 29] . The real 
space images shown in Fig. 3 is strikingly similar to the 
recent Bitter decoration of NbSe2 in the PE region [34) . 
which lends support to the idea that a phase transition 
is not necessary for the PE to occur. The amorphous 
state at b p is also consistent with the neutron scatte ring 
experiment 22] and the STM imaging of the vortices |5q 
in the PE region. 



D. Defect density rid(b) 

Figure 4(a) shows the defect fraction fd{b), defined as 
the ratio of number of defects to number of vortices, for 
A = 0.02. As mentioned earlier, reducing / to zero from 
I 3> I c generates a configuration which is stable and 
does not show hysteresis 47] . The fd(I = 0) can thus be 
considered as the equilibrium defects which determines 
the behaviour of stable I c {b). The /<j is typically greater 
than 0.35 for b < 0.1. Similar value for fd is reached at 
the onset of PE. The intermediate polycrystalline state 



(a) 

A = 0.02, N = 4096 
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FIG. 4: (a) The defect fraction f d (b) at I = 0. (b) The defect 
density n d (b) at I = 0. The n discl = \n d>7 - n<i,b\ is the 
approximate density of free disclinations, where n d ,7 and n dl 5 
are the density of vortices with coordination number 7 and 5, 
respectively. Inset: the n d (b) plot for A = 0.04. The thick 
dashed line is a fit to n d (b) = n d o exp(ab k ). 



is sandwiched between 0.3 < b < 0.7, with fd less than 
0.2. 

Interestingly, if the defect density rid (6) rather than 
fd(b) is plotted, a good fit can be obtained for b > 0.4 
using the expression 



rid (6) = n do exp(o!6 fc 



(3) 



This is shown in Fig.4(b) for A = 0.02 and A = 0.04 
(inset). The exponent k ks 2.6 and 2.4 for A = 0.02 and 
0.04, respectively. Surprisingly, both ndib) and I c {b) can 
be fit using the same value of k. This allows one to write 



Ic(b) 



Id 



1 



[n d {b)Y (i - 



(4) 



where r rj 2.5 and 2.7 for A = 0.02 and 0.04, respectively. 
Figure 4(b) also shows the density of free disclinations 
which can be approximately defined as ndiscl — \ n d.7 — 
rid,5 |j where n<j,7 and rid, 5 are the density of vortices with 
coordination number 7 and 5, respectively. The ndiscl is 
for b < bp^on, and rises rapidly in the PE region. 

From the analysis of the n d (b) and the I c (b) data, it 
is possible to conclude that I c {b) oc [rid(6)]~ r , where the 
effect of the pinning potential is implicit through the de- 
fect density ndib) and the exponent r. This empirical 
relation between n<j(6) and I c {b) captures the essential 
physics governing the behaviour of the critical current in 
a disordered vortex system with point impurities. Physi- 
cally, the increase in nd with b enables plastic shearing to 
be initiated at lower currents, which explains the mono- 
tonic decrease in I c below the PE. But this also leads to 
a paradoxical situation: since the nd(b) continues to rise 
rapidly even in the PE region, the enhanced plasticity 
of the system should cause the I c to decrease continu- 
ously till B C 2, contrary to the observation of PE above 
bp.on- The problem is compounded further since the real 
space images do not suggest any noticeable configura- 
tional change across b pon which can be associated with 
sudden increase in pinning due to quenched impurities. 

The above paradox can be resolved by conjecturing 
that the vortex system undergoes jamming at the onset 
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of the PE due to increase in dislocation density. The jam- 
ming envisaged here is similar to the jamming in gran- 
ular materials where the system develops a finite yield 
stress at a threshold packing fraction 56] . The jamming 
have been also studied in dislocation network formed in a 
plastically deformed material I57| which also shows com- 
plex spatio-temporal dynamics |58j|. The jamming in the 
vortex system could occur due to disordered packing of 
vortices which imposes kinematic constraint for the mo- 
tion of vortices. From Fig. 4(a), the defect fraction fa 
becomes 0.30 at b p on , which could possibly be the thresh- 
old value for the jamming to occur in the vortex system. 
As mentioned earlier, the dynamics close to the I c is gov- 
erned by few active channels. The effect of disordered 
packing of vortices is to increase the effective energy bar- 
rier AEt for the transverse mobility of the channel, thus 
requiring larger current to initiate channel dynamics. Fu- 
ture work should be able to obtain the -B-dependence of 
AEt, which could provide evidence for jamming in the 
PE region. 




~ <b) ' 










/ b = 0.60 
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FIG. 5: (a) The I c (b) and rid (ft) (inset) without the field de- 
pendence of the length scales A and £. The A = 0.02 and 
N v = 4096. (b) The 7 C and rid as a function of N v for b — 0.6 
and A = 0.02 



E. Field independent A and £ 

The I c (b) and the rid{b) obtained using X(b) = Ao and 
£(6) = £o are shown in Fig. 5(a). The Ao and £o are the 
values at B = 0. For b < 0.6, the I c (b) does not show any 
appreciable deviation from that of Fig. 1. At b rs 0.6, the 
defect density rid becomes zero, and the system remains 
topologically ordered till B C 2. As a consequence, the PE 
is absent at high fields and is consistent with the idea 
that the PE cannot occur in rigid vortex lattice. 



F. Finite size effect 

The finite size effect on rid and I c was determined by 
increasing N v for b — 0.6. The total number of defects 
shows a minimum around b = 0.60, or equivalently, the 
average size of the ordered region (domain) is maximum 
for this field. Hence, any relevant length scale determin- 
ing the system behaviour would be largest for b = 0.6. 
Figure 5(b) shows I c and rid as a function of the system 



size (oc N v ). The maximum N v which could be simu- 
lated was 6400. For N v > 4000, the I c tends to become 
independent of the system size though rid continues to in- 
crease, suggesting that the two quantities are determined 
by different length scales. The finite size calculation also 
shows that the system size used is sufficient to obtain 
macroscopic behaviour of the I c (b), and demonstrates the 
necessity of going to larger system size for doing realistic 
simulation of the vortex dynamics. 



IV. DISCUSSIONS 



The results presented above clearly demonstrates that 
a phase transition is not essential for the PE to occur in 
superconductors. This rules out any generic explanation 
based on order-disorder (or Bragg glass to vortex glass) 
transition, and resolves the difficulty in interpreting the 
PE in thin amorphous films. Also, the simulation results 
cannot be explained within LO theory. The R c in the ab- 
sence of dislocations in 2D was shown to be ~ 4 — 5an|53||. 
which is much smaller than the length scale determining 
the behaviour of I c . Further, the amorphous regime also 
exists for b < 0.2 and occurs due to single vortex pinning. 
If the amorphous state in the PE regime is attributed to 
the single vortex pinning, one expects similar dynamical 
behaviour in the PE regime and in low fields. Contrar- 
ily, the anomalous dynamics is observed only in the PE 
regime, thus indicating that one need to go beyond the 
pinning concept to describe the vortex state in the PE 
regime. 

As discussed in Section 3.4, an alternate scenario based 
on jamming can be invoked to explain the increase in I c 
with increasing rid- One of the important feature of a 
jammed state is the athermal relaxation under applied 
stress. For example, logarithmic relaxation is observed 
as a function of the number of "taps" applied to a a 
granular heap|59|. Evidence for similar behaviour exists 
in vortex system in the PE region. Henderson et aL|l4j| 
showed that the voltage response with bi-directional cur- 
rent grows with the number of pulses before saturating, 
which is an evidence for athermal behaviour. In an inter- 
esting recent experiment Ravikumar et aZ.[6(j measured 
the magnetic relaxation in VaSi in the PE region. The 
authors found that the magnetic moment M as a function 
of serial number of measurement n (rather than the real 
time) follows logarithmic behaviour. If the measurement 
n is assumed to be the number of "taps" applied to the 
system |6lJ, the athermal relaxation observed is similar 
to that observed in Ref . | 59| , indicating a common under- 
lying mechanism in the two systems. More experiments 
and simulations are necessary to probe the jamming sce- 
nario in the vortex system in the PE region. 
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V. CONCLUSIONS 

The paper describes the results from detailed simula- 
tion of the magnetic field dependence of I c for a super- 
conductor at T = 0. The I C {B) shows a large increase 
near the upper critical field B c2 in agreement with the 
experimental observation of the peak effect phenomenon. 
There is no evidence of a phase transition in the peak ef- 
fect region. Instead, a continuous transformation of a 
less disordered polycrystalline state into an amorphous 
state is observed implying that an ordered state is not 
necessary for the PE to occur. A jamming scenario is 
presented which could explain the PE and can account 



for the anomalous dynamics in the peak effect region. 
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